pub fn new_birthday_probability(n: u32) -> f64 {
  const DAY: f64 = 365.;
  let mut p = 1.;
  for i in 1..n {
    p *= (DAY - i as f64) / DAY;
  }
  1. - p
}
